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Inaugural Message from Editor-in-Chief on the Behalf of Editorial Board 


Dr. Nitish Puri, Dr. Ashwani Jain, Dr. Anil Kumar Sachan, Dr. P. K. Agarwal, Dr. S. K. Singh, Dr. K. K. Singh, Dr. Athar 


Hussain, Dr. Priyaranjan Pal, Dr. Pijush Samui, Dr. Sunita Kumari, Dr. Hemant Chore, Dr. Georgios Nikitas, Dr. Piyush 
Mohanty, Dr. Vishal Puri 


On behalf of the Editorial Board, it is our privilege to present the inaugural issue of the Annals of Earth, 
Environment, and Infrastructure Engineering, caringly shortened as AEEIE. We are enthusiastic to 
bring a change in the world of open-access publishing and aim to be a role model for open-access 
journals in India. We have launched this journal to provide good opportunities to students, researchers, 
and academicians to get their work reviewed, improved and published at low article processing charges 
without sacrificing the quality of the content and publishing ethics. The seasoned academicians from 
reputed institutes and research institutes and laboratories have been added to the editorial board. 


Annals of Earth, Environment, and Infrastructure Engineering is an open-access electronic journal. 
Original peer-reviewed articles are published in the journal in different categories like research papers, 
review papers, technical notes, case studies, and conference proceedings. The journal covers a broad 
spectrum of topics in Civil Engineering such as Geotechnical Engineering, Rock Mechanics, Soil 
Dynamics, Engineering Geology, Geophysics, Earthquake Engineering, Environmental Engineering, 
Transportation Engineering, Structural Engineering, Structural Dynamics, Water Resource 
Engineering, Irrigation Engineering, and other related inter-disciplinary fields. 


The inaugural issue has 5 papers, details of which are as follows: 
Assessment of Seismic Hazard for the District of Jalandhar (India) in view of Smart City Mission 
A Review of Failure Mechanism of Pile-supported Structures due to Earthquake Induced Soil 
Liquefaction 
Utilization of Municipal Solid Waste (MSW) Materials in Various Civil Engineering Works 
Leachate Recirculation Using Tire Scraps — An Economical Solution for Solid Waste Management 
Preliminary Assessment of the Seismic Hazard for the District of Pathankot (Punjab), India 


Launching this journal would not have been possible without the valuable support and suggestions from 
the Editorial Board members. We would like to take this opportunity to acknowledge all the authors 
who contributed to the inaugural issue of the journal. We are gratified by the esteemed members of our 
editorial board and anonymous reviewers for their time and effort in reviewing the manuscripts and 
providing valuable suggestions. Finally, we cordially welcome readers, researchers, and practitioners 
to submit manuscripts to this journal. Your suggestions and criticism would be welcomed at all the 
stages. 


Thank you all. We look forward to your positive response and encouragement. 
Dr. Nitish Puri 


Editor-in-Chief 
AEEIE (India) 
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ABSTRACT 


Probabilistic seismic hazard analysis (PSHA) has been carried out for the district of Jalandhar. The earthquake 
data has been collected from different seismological agencies, e.g. NDMA, IMD, ISC-UK and USGS, and a 
comprehensive earthquake catalogue has been developed. Several ground motion prediction equations have been 
reviewed to select a GMPE appropriate for carrying out PSHA. Earthquake hazard parameters have been 
estimated for PGA and spectral acceleration. The seismic hazard maps of district have been developed for 
different return periods at 10%, 2% and 1% probability of exceedance. The results show that the district can 
experience strong shaking from earthquakes originating in north-west Himalayas. Nonlinear site response analysis 
has also been carried out for limited number of sites in order to investigate the effect of local site conditions on 
earthquake ground motions. It has been observed that the soils of district of Jalandhar are capable of amplifying 
earthquake ground motions. Therefore, structures must be designed using the hazard parameters determined 
considering the local tectonic setup of the region as well as the local soil conditions. The developed hazard maps 
would help engineers and architects involved in planning and design of earthquake resistant structures and 
retrofitting works. 


Keywords: PSHA, Maximum magnitude potential, Peak ground acceleration, Spectral acceleration, Jalandhar 


INTRODUCTION of Delhi and Chandigarh city in the Union Territory 


Chandigarh. The Smart City Mission is a centrally 


Smart cities keep track of the condition of critical sponsored scheme financed by the Government of 


infrastructure facilities, which include water supply India with funds up to 2480 billion in a duration of 


network, rain water and waste water drainage five years. Retrofitting is the primary strategic 


systems, fuel supply and storage systems, electricity component of area-based development in the Smart 


supply units, communication systems, City Mission. Structures that need to be retrofitted 


transportation systems and building services, to in a city would be identified to build intensive 


better optimize resources, plan maintenance infrastructure facilities with smart applications 


activities, arrange immediate rescue in case of both packed in a Smart City. Planning of an efficient 


natural and man-made hazards, and monitor security retrofitting in a city requires good knowledge of 


aspects while maximizing services to citizens. In various expected hazards, e.g. flood hazard, wind 


India, a number of smart cities have been planned in hazard and earthquake hazard. 
various states, e.g. Karnal and Faridabad in the State 
of Haryana, Ludhiana, Jalandhar and Amritsar in the The present study focuses on_ the 


State of Punjab, New Delhi in the Union Territory assessment of earthquake hazard for the district of 


@ Corresponding Author’s Affiliation: 

AECOM India Private Limited 
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Jalandhar in State of Punjab (India) in view of 
planning for the proposed smart city ‘Jalandhar’. 


The earthquake zoning map of the India 
reported in IS 1893-Part1(2016), offers a broad 
zoning of the country for seismic hazard, and it gets 
updated only when new earthquakes occur. Hence, 
it is necessary to obtain reliable estimates of seismic 
hazard parameters calculated by carrying out state- 
of-the-art analysis on possible earthquake scenarios. 
The results of such analyses are formulated 
generally in terms of peak ground acceleration 
(PGA) and spectral acceleration (Sa) at various 
return periods. PGA is used to quantify earthquake 
ground motions, horizontal forces and shear stresses 
and strains in the earthquake resistant design 
procedures. Whereas, Sa is a crude representation of 
the response of a structure represented by the 
maximum acceleration experienced by single- 
degree-of-freedom system undergoing damped 
vibrations. 


The local soil conditions in top 20 to 30 m 
of soil profile play an important role in defining the 
characteristics of earthquake ground motions 
reaching at the surface of the earth. Site response 
analysis is often carried out to assess the effect of 
local soil conditions on_ earthquake wave 
propagation. The results of ground response analysis 
are expressed in terms of amplification factors for 
PGA and Sa. The values of PGA and Sa for rock 
sites, obtained from either deterministic or 
probabilistic seismic hazard assessments can be 
suitably modified using the amplification factors 
obtained from site response analysis. 


In the present study, probabilistic seismic 
hazard analysis and nonlinear site response analysis 
has been performed for the district of Jalandhar, and 
seismic hazard maps based on PGA and Sa have 
been prepared at various return periods. The 
developed hazard maps along with the suggested 
amplification factors would help engineers and 
architects involved in planning and design of 
earthquake resistant structures and retrofitting 
works in the District. 


STUDY AREA AND TECTONIC MAP 
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The entire Jalandhar Division was awarded to India 
when Punjab was partitioned. Jalandhar District, 
with population of about 2.19 million is the 4" most 
populated district in the state of Punjab. The district 
of Jalandhar is the 9" biggest district in terms of area 
in the state of Punjab with an approximate 
geographical area of about 2624 km’. The 
population density in the district 1s approximately 
836 people per km”. 


The region has experienced some mild 
seismic activity in the past in the range of M 4.0-5.0. 
On 14" March 2010, the region experienced a mild 
earthquake of magnitude of 4.5, that occurred in 
northern Punjab along the Punjab-Himachal Pradesh 
border. It was felt over a wide area due to its shallow 
depth near Hiranagar (Punjab). In August 2013, 
another mild intensity earthquake of magnitude 4.7 
on the Richter’s scale was felt at many places in 
Punjab. The epicenter of the quake was in the 
Hoshiarpur-Himachal Pradesh border region as per 
India Metrological Department (IMD). However, 
some disastrous earthquakes have also occurred near 
the region in the last 200 years, which are 1934 
Bihar-Nepal earthquake (Mw 8.4), 1950 Assam 
earthquake (Mw 8.7), 1905 Kangra earthquake (Mw 
7.8), 1991 Uttarkashi (Mw 6.8), 1993 Killari 
earthquake (Mw 6.2), 1999 Chamoli earthquake 
(Mw 6.8,), 2005 Kashmir earthquake (Mw 7.6) and 
2015 Gorkha earthquake (Mw 7.8). The havoc 
produced by these earthquakes has been a wakeup 
call for the government to take suitable mitigation 
measures. 


An area of 300 km radius around Jalandhar 
district has been selected as the seismic study region 
(28-34 Degree N to 72-78 Degree E). A tectonic 
map developed for the seismic study region by Puri 
and Jain (2018) has been adopted for the study (Fig. 
1). The tectonic map was developed using SEISAT 
(Dasgupta et al., 2000). The potential earthquake 
hazard in the study area is governed by several 
tectonic features, which include Himalayan Frontal 
Thrust on the north and north-eastern side, local 
lineaments on the south and Sargodha-Lahore-Delhi 
ridge on the north-western side. 


ESTIMATION OF GUTENBERG-RICHTER 
PARAMETERS 


The development of a comprehensive catalogue of 
seismic events is a first step for any seismic hazard 
assessment program. A_ reliable estimate of 
earthquake hazard in a region strongly depends on 
the accuracy and adequacy of the information on 
earthquake events. Information on damaging 
earthquakes that have occurred during pre- 
instrumental period has been collected from the 
catalogue developed by National Disaster 
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Management Authority of India (NDMA). The 
earthquake data for the instrumental period have 
been collected from various international and 
national earthquake monitoring agencies, e.g. 
United States Geological Survey (USGS), India 
Meteorological Department (IMD) and National 
Disaster Management Agency (NDMA\), India 
(NDMA, 2011), and International Seismological 
Centre, UK (ISC-UK). An earthquake catalogue has 
been developed, for 1800 AD to 2017 AD, by 
combining the pre-instrumental and instrumental 
data. 


Fault 
Minor Lineament 
Thrust 
Ridge 
Fault Involving Basement and Cover 
Major Lineament 
Subsurface Fault 
Suture 

- Study Area 
6-6.9 Mw 
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Jalandhar 





Fig. 1. Tectonic setup of the seismic study region (after Puri and Jain, 2018) 


Table 1. Completeness analysis of catalogue 


Magnitude class 


(Mw) 
Period 
4.0-4.9 1962-2015 
5.0-5.9 1926-2015 
6.0-6.9 1901-1975 
7.0-7.9 1905-1999 
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CUVI method 


Completeness 
Stepp method 
Interval (Years) Period Interval (Years) 
53 1963-2015 a2 
89 1925-2015 90 
74 1900-1970 70 
94 1904-1999 95 


The events with earthquake magnitudes 
other than My have been converted using the 
available empirical relationships between various 
earthquake magnitude scales and moment 
magnitude (M,,). Pre-instrumental data is usually 
available on MMI scale and has been changed to 
moment magnitude by Gutenberg-Richter equation, 
which is as follows: 


My = 2/3 x MMI + 1 (1) 


The correlations developed by Scordilis (2008) 
between mp - My and M, - My have been used. 


The following correlation derived by Kolathayar et 
al. (2012) has been considered for the conversion of 
local magnitude to moment magnitude. 


My = 0.815 x M, + 0.767, for3.3<Mri< (2) 
7.0 


The following correlation developed by Yenier et al. 
(2008) between Mg - My has been used. 


My = 0.764 x Ma + 1.379, for3.7<MaS (3) 
6.0 

An epicentral map has been developed using the 
prepared catalogue (Fig. 2). The catalogue has been 
checked for completeness using CUVI (Tinti and 
Mulargia, 1985) and Stepp method (Stepp, 1972) 
and found to be complete for a sufficient period of 
time. The completion periods have been reported in 
the Table 1. The compiled catalogue contains 1256 
earthquake events of moment magnitude (My) > 4 
that have occurred up to December 2017 in the study 
region. 


For identifying potential seismogenic 
sources, the epicenters from the compiled catalogue 
(Fig. 2) have been plotted over tectonic map (Fig. 
1). A total of 10 tectonic features associated with 
earthquakes of My = 4.0 have been identified as 
active seismogenic sources, out of these, 8 sources 
are potential seismogenic sources as given in table 
below. 

The earthquakes of My => 5 have been 
considered for Mmax estimation. The Mmax values for 
the major seismogenic sources in the study area have 
been estimated by adding an increment 0.5 to the 
respective values of Mobs (Gupta, 2002). The 
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calculated Mmax values for study region are between 
moment magnitude of 5.5 to 8.5 (Table 3). 


The Gutenberg-Richter 
parameters a and b are the prime input parameters in 
PSHA and hence need to be calculated with 
precision. For this, the study region was divided into 


seismicity 


three sub-regions considering each sub-region as an 
area source of earthquakes. The seismicity 
parameters were computed for each area source 
through linear least squares regression method 
following an exponential distribution of magnitude 
taking into account the complete part of the 
catalogue for all the magnitude classes. The 
exponential distribution is given in equation (4) 
below: 


Am = 102->!w = exp(a — BM,,) (4) 


where, Am = mean yearly rate of exceedance, a = 
coefficient such that a power of 10 gives the mean 
annual number of earthquakes of magnitude > 0, a = 
2.303a, b = coefficient that describes the possibility 
of large and small earthquakes and B = 2.303b. 


The reciprocal of Am for a given magnitude 
is called return period (Tr) of an earthquake 
exceeding that magnitude and is very important for 
earthquake resistant design. The seismicity of all 
area sources considered in the study has been 
calculated and reported in Table 4. The value of 
return period calculated for different area sources 
demonstrates the capability of tectonic sources in 
Himalayan Thrust System to generate frequent large 
earthquakes. 


Assesement of Seismic Hazard 


The global ground motion prediction equation 
developed by Abrahamson et al. (2013) has been 
used for the estimation of earthquake hazard. PSHA 
has been done considering three sub-regions, viz. 
Himalayan Thrust System, Sargodha-Lahore-Delhi 
Ridge and Aravalli-Delhi Fold Belt as area sources 
of earthquakes and based on the developed 
catalogue, average focal depths have been taken as 
15 km, 17 km and 10 km for the selected three sub- 
regions, respectively. 
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Fig. 2. Epicenter Map of Study Area 
Table 2. Maximum Observed Magnitudes for Various Sources 
S. No. Seismogenic Source Mobs 

I Jwala Mukhi Thrust (JMT) 5.5 
23 Kaurik Fault System (KFS) 6.8 
3. Lineament System of SLDR (LSLDR) 6 
4. Main Boundary Thrust (MBT) 8 
> Main Frontal Thrust (MFT) 5 
6. Main Crustral Thrust (MCT) ips 
a Sundar Nagar Fault (SNF) 7 
8. Sargodha Lahore Delhi Ridge (SLDR) 6.5 


Table 3. Minax values for Various Sources 


Seismogenic Source Total Length of Fault (in km) Minax 
Sargodha Lahore Delhi Ridge (SLDR) Area Source 7 
Main Boundary Thrust (MBT) 825 8.5 
Main Frontal Thrust (MFT) 32 5.5 
Jwala Mukhi Thrust (JMT) 290 6.0 
Main Crustral Thrust (MCT) 769 7.8 
Sundar Nagar Fault (SNF) 101 7.5 
Kaurik Fault System (KFS) 137 (Pe 
Lineament System of SLDR (LSLDR) 97 6.5 
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Table 4. Seismicity parameters for different area sources 


Area source b 
Himalayan Thrust System 0.75 
Aravalli-Delhi Fold Belt 0.69 


Sargodha-Lahore-Delhi Ridge 0.85 


The hazard has been calculated for 10%, 
2% and 1% probability of exceedance in time frame 
as 50 years as per recommendations of Eurocode 8: 
BS-EN 1998-1 (2005). For ordinary structures, the 
seismic hazard map for 10% probability of 
exceedance in 50 years is recommended. However, 
for important facilities like Nuclear Power Plants 
and other megastructures, the maps for 2% and 1% 
probability are used. The PSHA software R-CRISIS 
v. 18.2 (Ordaz and Salgado-Galvez, 2017) has been 
used for the purpose. R-CRISIS allows the 
calculation of results for exceedance probability 
plots, set of stochastic events etc. The parameters 
such as a, b, Mmin, Mmax, Am and attenuation models 
are the input parameters, and PGA and PSA are the 
outputs. 


The hazard maps have been prepared for 
return periods of 475 years, 2475 years and 4975 
years at 10%, 2% and 1% _ probabilities of 
exceedance in 50 years respectively (Fig. 3 to Fig. 
4). The estimated peak ground acceleration ranges 
from 0.08¢ - 0.13, 0.15g - 0.26g and 0.23¢ - 0.34¢g 
at 10%, 2% and 1% probabilities of exceedance 
respectively. It has been observed that the IS code 
underestimates the hazard at 2% and 1% probability 
of exceedance. The response spectra have been 
constructed for various return periods corresponding 
to maximum observed PGA and shown in Fig. 5 
along with those specified in IS 1893 Part-1 for 
comparison. Hazard maps based on_ spectral 
acceleration (g) at 0.1 sec, 0.2 sec, 1.0 sec and 2.0 
sec at 10%, 2% and 1% probabilities of exceedance 
have been also developed (Fig. 6). It has been 
observed that the IS code underestimates the 
spectral acceleration at higher structural periods for 
all return periods. Hence, high rise buildings in these 
areas constructed as per current Indian Standards 
could receive damage during earthquakes. However, 
a detailed study on the effect of local site conditions 
would provide a deeper insight into the problem at 
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a Magnitude (Mw) class 


3.8 4.0-8.0 
2.61 4.0-7.0 
3.27 4.0-6.5 


hand. 
1D Site Response Analysis 


Site response analysis has been carried out for 5 sites 
(Fig. 7) drilled up to refusal in District of Jalandhar. 
For this, one dimensional nonlinear model has been 
adopted to estimate the response using DEEPSOIL 
software (Hashash et al., 2016). The methodology 
has been adopted from Puri and Jain (2018). Based 
on geotechnical data collected, sites are classified as 
class D sites by calculating average SPT N-value of 
the profile as per the recommendations of NEHRP 
(FEMA 368, 2000). 


The stiffness and damping of soil layer 
play fundamental role in_ estimating wave 
amplification parameters in seismic microzonation 
studies. The analysis requires characterization of the 
stiffness of an element of soil considering low strain 
Shear modulus (Gmax), variation of modulus ratio 
(G/Gmax) and Damping (D) with cyclic strain 
amplitude (vy), and other parameters. Due to 
unavailability of in-situ Shear wave velocity (Vs) or 
Shear Modulus (Gmax) recordings, a correlation 
between Gmax - N developed by Ohsaki and Iwasaki 
(1973) for sandy soils has been used and shear 
modulus degradation (G/Gmax-y) and damping ratio 
(D-y) curves developed by Darendeli (2001) have 
been used. The Darendeli curves require the value 
of coefficient of earth pressure at rest (K.) for each 
layer in order to fit G/Gmax - y curves and D-y curves 
for each layer. The K, values have been calculated 
using the following equation: 


K, =1-singd (5) 
where = angle of internal friction. The values 
have been estimated from SPT N-values using the 
correlation developed by Puri et al. (2018). The 
typical input parameters for site 4 have been 
reported in Table 5. 
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Fig. 3. Seismic Hazard Map of Jalandhar at 10% Probability of Exceedance in 50 years 
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Fig. 4. Seismic Hazard Map of Jalandhar at (a) 2% Probability of Exceedance in 50 years and (b) 1% Probability of Exceedance in 50 years 
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Fig. 5. Comparison of PSHA based response spectra with response spectra specified in IS 1893-Part 1:2016 
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Suitable acceleration time histories can be 
selected based on magnitude of controlling 
earthquake, PGA value, source to site distance and 
site class. PGA values obtained from PSHA for rock 
sites corresponding to a return period of 475 year 
have been used for the selection of input motions at 
each site. Acceleration time history of the 1991 
Uttarkashi earthquake M,, 6.8 (focal depth = 10 km) 
recorded at the Ghansiali station with PGA = 0.118¢g 
is used for the analysis (Fig. 8). 
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Fig.7. Location of boreholes used for site response analysis 


Table 5. Typical input parameters for Site no. 4 


Depth Soil PI Unit SPT 

Type weight 'N' 

(kN/m?) 

1.5 Sand NP 15.37 7 
3 Sand NP 15.65 10 
4.5 Sand NP 15.84 12 
6 Sand NP 17.16 26 
Too Sand NP 17.72 32 
9 Sand NP 18.20 37 
10.5 Sand NP 21.90 61 2 
12 Sand NP 23.56 73 = 
13.5 Sand NP 25.77 89 g 
15 Sand NP 26.19 92 
16.5 Sand NP 27.15 99 
18 Sand NP 20.38 50 
19.5 Sand NP 20.38 50 


*NP = Non-plastic; K, = Coefficient of earth pressure at rest 


Acceleration (g) 





0 5 10 15 20 25 30 35 40 45 
Time (seconds) 


Fig.8. Input acceleration time history of 1991 Uttarkashi 
earthquake recorded at Ghansiali station 


The results of the wave amplification analysis have 
been reported in Table 6. Fig. 9 shows a comparison 
of PGA (g) and shear strain (%) observed at 
different sites. Due to limited borehole data (eight 
boreholes), interpretation cannot be made for the 
amplification trend across the District. However, 
based on significant amplification observed at five 
sites, it can be concluded that the city may 
experience high ground accelerations. Moreover, 
high strains are observed for all the five sites and 
there is a possibility of substantial settlements 
during an earthquake. The amplification factors for 
the analyzed sites range from 1.65 to 2.11 with an 
average value of 1.9. The maximum PGA of 0.249¢ 
and minimum PGA of 0.195g are observed for 
Bulandpur and Wazir Singh Enclave sites, 
respectively. Based on the average amplification 
factor observed for PGA, the region can experience 


PGAs ranging from 0.152g to 0.247g. However, a 


Annals of Earth, Environment, and Infrastructure Engineering 


Volume 1, Issue 1 (November-December 2021), pp. 1-12 


N vs Depth Gmax Gmax vs Depth Ko 
(Mpa) 
0 55.51 0 0.51 
5 65.89 2 0.50 
‘4 136.30 4 0.43 
F 165.68 6 0.41 
189.90 0.38 
8 ie 
303.83 = 0.29 
10 = 10 
359.71 2 0.24 
v A 12 
433.37 0.19 
14 447.08 a 0.18 
16 478.99 16 0.16 
FP 7 252.03 18 0.33 
252.03 20 0.33 


20 


more detailed study specifically on site-response 
analysis is required for the region. The amplification 
factors for spectral acceleration (average of all sites) 
have been shown in Fig. 10. It has been observed 
that at all structural periods the spectral acceleration 
has been amplified. This shows that the estimated 
seismic scenario for the District of Jalandhar is 
worse than that proposed by the Indian Seismic 
Code. Hence, at higher structural periods, severe 
dynamic loading could be inflicted to structures 
designed using current structural design and 
of India. The 
amplification factors and the hazard maps presented 


construction bylaws average 


in this study can be used in design and retrofitting of 
structures in district. 


CONCLUSION 


The assessment seismic hazard has been carried out 
using probabilistic approach for estimating hazard 
parameters for rock sites for district of Jalandhar. 
The district is always under the threat from 
earthquakes due to its proximity to Himalayan 
frontal fault. It has been observed that the district is 
prone to high to severe earthquake ground motions 
ranging from 0.152g to 0.247g. Hence, a 
comprehensive estimation of earthquake hazard for 
the region is required. The results of this study 
would help and motivate engineers and investigators 
involved in seismic design and retrofitting of 
structures and planning for mitigation measures. 
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Table 6. Results of site response analysis 


Site No. Latitude (N) Longitude (E) 
1 31.38363 75.54201 
2 31.28898 75.59628 
3 31.34913 75.55048 
4 31.37636 75.5883 1 
5 3132522 75.63357 
PGA (g) 
0.05 0.1 0.15 0.2 0.25 


bo 


Depth (m) 


10 


—— Site 1 
—?@ Site 2 
—t— Site 3 
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14 


16 
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Kf 


Amplification Factor for Sa (g) 


0.01 0.1 


Surface PGA (g) 
0.205 
0.244 

0.23 
0.249 
0.195 


Depth (m) 


Period (sec) 


Amplification Factor 


1.74 
2.07 
1.95 
2.11 
1.65 


Strain (%) 
0.1 0.2 0.3 
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Fig.10. Observed amplification factors for spectral acceleration (Sa) 
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ABSTRACT 


The paper reviews the failure mechanism of pile-supported structures in liquefiable soil during earthquakes. Post- 
earthquake investigation of many pile-supported structures on liquefiable soil shows that despite the large factor 
of safety against bending due to lateral loads employed in their design, structures failed, e.g. failure of firehouse 
building at Kobe Port in 1995, Japan. It has also been observed that the middle piers of the bridges collapsed 
without significantly affecting the other parts. A number of experiments conducted by various investigators 
highlighted the importance of pile failure modes other than bending mechanism, especially, pile failure due to 
buckling. The initial codes of practice for pile design on liquefiable soils such as NEHRP 2000, Eurocode 8, 
JRA1996, and IS1893 focused only on bending failure, and failure due to buckling/dynamic amplification of pile 
due to the axial loading during liquefaction was neglected. 

Keywords: Pile failure, liquefiable soil, bending, buckling, shear, dynamic amplification 

INTRODUCTION movement) induces bending in the piles. Earthquake 
exerts an abrupt dynamic loading and causes 
Pile foundations are most widely used deep damage to structures both above and below the 
ground. It is very difficult to ascertain the failure 
mechanism of foundation of a structure without 


foundation especially in case of _ saturated 
cohesionless soil, which is_ susceptible to 


liquefaction during an earthquake. Due to the rise in deep excavation. Twenty years after the Niigata 


earthquake 1964 and the Kobe earthquake 1995, 
investigations have been done to find the failure 
pattern of the piles. Piles were excavated from the 
subsoil, photographs of boreholes were taken, and 
pile integrity tests were carried out. These studies 
indicated the location of the cracks and damage 


global population and limited land resources, 
infrastructures are often constructed over reclaimed 
lands e.g. the port Island in Kobe (Japan), Rajorhat 
and Saltlake in Kolkata (India). Such structures are 
prone to liquefaction induced damages if subjected 
to dynamic loading (earthquakes). Structures on 


loose to medium dense sands are often provided pile patterns for the piles. Some of the observations are 


foundations as the bearing soils are not stiff enough presented below: 


to support these structures. Failure of pile 


foundations was noticed in many cases and it was 1. The failure of foundations occurred in both 


observed that even after the provision of large factor 
of safety against bending failure, structures often 
failed. In bending failure mechanism, lateral loads 
due to inertia and lateral spreading (slope 


@ Corresponding Author’s Affiliation: 
COWI, Plot No 121, Phase I, Udyog Vihar, Sector 20, 
Gurugram, Haryana (India)- 122016 


laterally spreading grounds and level grounds. The 
failure pattern (tilting of buildings) is similar. This 
is shown in Fig 1. 
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Fig.1. Building Failures after Kobe Earthquake (after 

Bhattacharya et al, 2013) 


2. Cracks were observed near top and bottom 
boundaries of liquified and non liquified soil layers. 
Bhattacharya & Madabhushi (2008) observed the 
formation of plastic hinges at different depth of 
liquified and non liquified soil layers. 
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Fig. 2: SPT values describing variation of soil profile along the 
depth of pile 

3. It has been observed that large diameter piles 
performed well during the seismic activity. For 
example, the older jetties of the Kandla port 
supported on R.C.C. piles of diameter 0.5 m suffered 
cracking during the Bhuj earthquake 2001, whereas 
in the same area, the RCC piles having 1.0 m 
diameter on the newer jetties performed well. 


DIFFERENT TYPES OF FAILURE 
MECHANISMS FOR PILE SUPPORTED 
STRUCTURES IN LIQUEFIABLE SOILS 


1. Bending failure mechanism: Bending failure 
generally occurs due to transverse or lateral seismic 
load or due to inertia of the superstructure. As the 
earthquake excitation begins, soil near the pile starts 
losing strength and flows in lateral direction, which 
is also responsible for the bending failure of the pile. 
Effect of this lateral flow of soil will be more in case 
of sloping ground called lateral spreading of the 
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ground. Bending in the pile due to lateral spreading 
of ground is often regarded as the root cause of many 
bridge failure. 


2. Buckling failure mechanism: As the process of 
liquefaction progresses, confinement of pile by 
surrounding soil decreases to such a low value that 
the pile starts acting like a long slender column. 
Now, axial load coming from the superstructure will 
cause buckling and pile will not be able to support 
the superstructure load. Magnitude of buckling will 
depend upon the stiffness of the pile, depth of the 
liquefiable soil and magnitude of ground shaking 
during earthquakes. 


3. Shear failure: This type of failure occurs due to 
lateral load such as kinematic or inertia load or 
combination of both. This type of failure is generally 
observed in hollow circular concrete piles having 
low shear strength capacity. 


4. Settlement failure mechanism: Due to combined 
action of buckling and bending it may be possible 
that structure will settle vertical or it may tilt in one 
direction. 


5. Dynamic failure mechanism: At the central pier 
of the river bridge, the depth of water increases in 
comparison to the water level near abutment and 
depth of embedment decreases and hence, the 
unsupported length of pier increases. This elongates 
the natural period for the central bridge pier as 
compared to the neighboring piers and abutment. 
Many bridge failures are attributed to collapse of 
central piers. The collapse of Showa Bridge (1964 
Niigata earthquake) is considered to demonstrate 
this mechanism. Table 1 summarizes the list of 
various bridge failures due to liquefaction. Table 2 
summarizes the performance review of various 
bridges situated in earthquake impacted locations. 





Fig. 3: Various Failure Mechanism a) Kandla Port Building 
with pile foundation; b) shear failure; c) Bending failure; d) 
Buckling; e) Dynamic amplification (Bhattacharya et al., 2008) 
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Table 1. List of Bridge Failures due to Soil Liquefaction (Bhattacharya et al., 2017) 


Earthquake Event Bridge Name Damage Type 
Tohoku (2011) Rokko Bridge Mid span collapsed 
Wenchuan (2008) Miaoziping Bridge One of the approaching spans collapsed 
Gaoyuancun Bridge Mid span fell off the piers due to liquefaction 
Costa Rica (1991) Rio Viscaya Bridge One internal supporting pier was missing 
Philippines (1990) Magsaysay Bridge Piers settled and failed, and the bridge fell into the river 
Tangshan (1976) Zhuacun Bridge The girders of the mid span collapsed 
Shahe Bridge Bridge girder support collapsed 
Philippines (1976) Quirino Bridge Mid span collapsed for the truss bridge 
Haicheng (1975) Panshan Bridge One of the mid piers sank, causing collapse 
Niigata (1964) Showa Bridge Few mid span collapsed 


Table 2: Assessment of Performance of Various Pile Supported Bridges during Earthquakes (Source: Bhattacharya et al., 2004) 


Case History and Reference Pile section/ type Performance Case History and Reference Pile section/ type Performance 
10 Storey Hokuriku Building, Hamada (1992a, b) 0.4 m dia RCC Good N.H.K Building, Hamada (1992a, b) 0.35 m dia RCC Poor 
hollow 
Landing Bridge, Berrill et al. (2001) 0.4 m square PSC Good NFCH Building, Hamada (1992a, b) 0.35 m dia RCC Poor 
14 Storey Building, Tokimatsu et al. (1996) 2.5 m dia RCC Good Yachiyo Bridge, Hamada (1992a, b) 0.3 m dia RCC Poor 
Hanshin Expressway pier, Ishihara (1997) 1.5 m dia RCC Good Gaiko Ware House, Hamada (1992a, b) 0.6 m dia PSC Poor 
hollow 
LPG Tank 101, Ishihara (1997) 1.1 m dia RCC 4 Storey fire house, Tokimatsu ef al. (1996) 0.4 m dia PSC Poor 
Kobe Shimim hospital, Soga (1997) 0.66 m dia steel tube Good Showa Bridge, Hamada (1992a, b) 0.6 m dia steel tube Poor 
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CASE STUDY OF DIFFERENT STRUCTURES 
DURING VARIOUS EARTHQUAKE 


The case studies reported by Dash et al. (2012) have 
been reviewed and presented below: 


1. Failure of Customs Office Tower at Kandla 
(India): The Customs office tower at Kandla port is 
a six-floor building, 22 m in height and very close to 
the waterfront (Figure 4). This pile-supported tower 
tilted about 30 cm at its top during the Bhuj 
earthquake 2001. The building was founded on 32 
short cast-in-situ concrete piles. Each pile was 0.4 m 
in diameter and 18m long. The piles were passing 
through clayey crust of 10 m depth and then 
terminated in a sandy soil layer. At this site, very 
little damage occurred in the super-structure, but 
considerable damage happened in the foundation. 
The analysis was carried out by applying the 
superstructure load and the lateral spreading 
simultaneously to investigate the combined 
response of lateral spreading and settlement. 


0.3m 1m 






«<< Before canhquake 
After earthquake 





0.3m 54 
' Boundary wall 


Sapper 8 Oe eee nines 
Brn tree tens 


Fig 4: Plausible settlhement mechanism of failure showing the 
tilting the Tower, assuming there is no structural failure of piles. 


2. Failure of firehouse building in Kobe Port 
(Japan): Four-storey firehouse building at the foot 
of the Kobe O-Hashi Bridge on Port Island moved 
and tilted towards the sea during Kobe earthquake 
(1995). Bhattacharya (2006) studied this case 
history to understand the mechanism of its pile 
failure. The piles were axially loaded heavily (412 
kN each) and were laterally unsupported for 16.4m 
of their length due to liquefied soil layer (Figure 5). 
The factor of safety against buckling for the piles 
was found to be 0.81. This suggests that the piles 
could probably have failed due to buckling at full 
liquefaction of the ground during the earthquake. 


Annals of Earth, Environment, and Infrastructure Engineering 


Volume 1, Issue 1 (November-December 2021), pp. 13-17 


Inclinaion 
building - + 
Quay 
G.L.0 
= = 
Inclination of pile 5 Fill 
=about 1/10 
Horizontal crack = 
Pile connections 


= Sandy sit === 





50 Tag Te hU el ee 


* ou he Ea eZe 


Fig. 5: Schematic of building’s foundation failure during 1995 
Kobe earthquake 


3. Failure of NHK Building in Niigata City (Japan): 
Four-storey R.C frame building supported over RC 
piles leaned during the Niigata earthquakes (1964). 
Widespread liquefaction and lateral spreading of the 
soil were observed near the building site due to 
which bending and dynamic amplification occurred. 


CONCLUSION 


The failure mechanism of pile foundation in soil, 
which is susceptible to liquefaction, is difficult to 
predict just based on surface observation. Engineers 
need to carry out detailed calculations. Similar type 
of surface observation is possible for different pile 
failure mechanisms. For example, bending, 
buckling, shear, uneven settlement and dynamic 
amplification of the pile group may lead to tilting or 
collapse of the superstructure. All the three 
structures studied in this paper were constructed on 
liquefied soils and titled after the occurrence of a 
damaging earthquake, however the governing 
mechanism of their pile failure were very different. 
It is, however, important to recognize the failure 
mechanisms, which will govern the design of new 
foundations or retrofitting of existing foundations. 
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ABSTRACT 


The effects of the addition of processed municipal solid wastes (MSW) materials on the engineering properties 
of M25 concrete and soft soils were investigated. For this, abandoned pile heads were collected and crushed into 
20 mm sized aggregates. PVC waste, leather waste, nylon fiber cuttings, and iron waste were collected from 
nearby industries and pulverized to the required particle size. The processed waste materials were used as a 
replacement to conventional concrete aggregates. Many improvements have been observed in terms of slump, 
weight, and compressive strength of resultant concrete. However, leather waste was proved to be unsuccessful in 
improving any engineering property of concrete. Materials like cement kiln dust (CKD), Flyash, and Rice husk 
ash (RHA) were utilized for stabilizing subgrade and foundation soils. Significant improvements in strength were 
observed for stabilized soil samples. However, only CKD was proved successful in improving the compaction 
characteristics of the soil. Based on the study, waste like recycled aggregates, nylon cuttings, and iron waste are 
recommended to be used as a replacement to conventional aggregates for concrete. Also, municipal solid wastes 
like CKD, Flyash, and RHA are recommended to be used as soil stabilizers but their possible effects on human 
health should be taken care of in advance. 


Keywords: Municipal solid waste, recycled aggregates, slump, compressive strength, soil stabilizers, unconfined 
compressive strength 


INTRODUCTION 


The world is passing through an era of urbanization 
and industrialization. Many infrastructure projects 
are being launched by the governments in every 
country. Simultaneously, numerous industrial units 
are being installed in almost every part of India as 
well. However, municipal authorities are facing 
huge difficulty in disposing of the industrial wastes/ 
construction debris being generated from these 
industries. Moreover, the existing dumpsites have 
also been running overloaded due to slow rate as 
well as poor technology available for recycling. This 
has made the issue of waste utilization a need of the 
hour. 
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Different types of waste materials are 
generated as a by-product of industrial operations 
and construction/ demolition activities. These 
materials can be utilized in various civil engineering 
works. In this study, different waste materials have 
been utilized to improve the engineering properties 
of cement concrete, pavement subgrades, and 
foundation soil. Coarse and fine aggregates in 
concrete were replaced with PVC waste, nylon fiber 
cuttings, iron waste, construction debris, and leather 
waste. Properties like weight, slump, workability, 
and compressive strength were determined for the 
resultant concrete blocks. On the other hand, 
materials like Rice husk ash (RHA), Fly ash, and 
Cement kiln dust (CKD) were used as additives to 
improve the properties of soft soils. Engineering 
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properties like shear strength and compressibility 
were assessed for the stabilized soil samples. 


In the last two decades, several 
investigators have utilized these MSW materials to 
augment the properties of cement, concrete, and soil. 
Sikalidisa et al. (2002) suggested a methodology to 
prepare mortar from municipal solid waste. An 
extensive piece of work is done on recycled 
aggregates by Claudio and Antonio (2011). It was 
observed that the compressive strength of concrete 
from recycled aggregates is marginally less than that 
of concrete having conventional aggregates. 
However, durability was observed to be the same. 
Nassar and Soroushian (2012) utilized milled glass 
as a partial replacement for cement in concrete. 
Significant improvement in strength and durability 
was observed and this improvement was attributed 
to the formation of secondary calcium silicate 
hydrate. An extensive review of RCA was published 
by McNeil and Kang (2013). They observed a small 
decrease in values of compressive strength, modulus 
of rupture, and modulus of elasticity when 
conventional aggregates were replaced by recycled 
aggregates. However, a significant improvement in 
tensile strength was observed. Also, _ the 
performance of the structure is marginally affected 
by the replacement of natural aggregates with 
recycled ones. It was recommended by them to use 
recycled aggregates in concrete structures. 


Stabilization is one of the oldest soil 
improvement techniques. It is being utilized since 
roman times. Dahale et al. (2012) published a 
review on soil stabilization. He concluded that 
utilization of waste materials in soil stabilization 
will not only solve the problem of waste disposal but 
will also provide new material alternatives. Moses 
and Saminu (2012) observed — significant 
improvements in compaction and __ strength 
characteristics of black cotton soil when it was 
treated with CKD. Recently Pandey et al. (2013) 
have utilized lime, jute fiber, and Fly ash 
individually and in combination to stabilize black 
cotton soil. They observed that a combination of the 
three waste materials was much more efficient in 
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stabilizing black cotton soils as compared to the 
waste materials alone. 


The purpose of this study was to see the 
effect of municipal solid wastes in improving the 
engineering properties of concrete and problematic 
soils. A better understanding of these characteristics 
will enhance the usage of these materials in 
construction works in places where they are in 
abundance. The study also focuses at the reduction 
of huge stockpiles of the various municipal solid 
wastes and their potential impact on human health. 


MATERIALS AND METHODS 
1. Concrete and Solid Waste 


a) Physical properties: The cement used was of 43 
grade and its physical properties are reported in 
Table 1. Abandoned concrete pile heads were 
collected. These pile heads were crushed into 20 mm 
size particles and used as aggregates in concrete. In 
this paper, these aggregates are referred to as 
recycled aggregates. Waste leather was collected 
from nearby leather industries and it was pulverized 
to fine aggregates-sized particles. Similarly, 
pulverized PVC waste and iron waste were used as 
a replacement for fine aggregates. For iron waste, 
fineness modulus, specific gravity, and density were 
observed to be 2.65, 4.5, and 1.97 g/cc respectively. 
The density of nylon fiber cuttings was observed to 
be 1.13 g/cc. Physical properties of different 
aggregates are reported in Table 2. All the properties 
were determined in the laboratory as per the 
guidelines of IS codes ([S:383, 1970; IS:516, 1959; 
IS:1199, 1959). 


b) Mix-design: Recycled aggregates were used to 
replace conventional aggregates in concrete. Three 
types of samples were prepared using 100%, 50%, 
and 25% recycled aggregates as replacements to 
conventional aggregates. PVC waste and leather 
wastes being non-earthen materials were used in 
small quantities. Fine aggregates were replaced by 
PVC waste and leather wastes in proportions of 5 
and 10%. Iron waste was used as a replacement to 
fine aggregates in the proportion of 20%. 
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Physical Properties 


Specific Gravity 


Water Absorption (%) 
Moisture content (%) 
Bulking (%) 

Fineness Modulus 


Grading Zone 


Tablel: Physical Properties of 43- Grade Cement 


Physical Properties 

Specific Gravity 

Fineness (%) 

Standard Consistency (%) 
Initial Setting Time (minutes) 
Final Setting Time (minutes) 
Soundness (mm) 


28 days Compressive Strength (N/mm7’) 


Observed Value 


3.07 
3.5 
34 


45 


Table 2: Physical Properties of Materials Used as Aggregates 


Sand Coarse Recycled Aggregates 

Aggregates 190% = 50% 25% 
2.67 3.01 2.66 2.83 2.92 
1.020 3.806 ye 4.33 4.21 
0.155 0.806 0.6 0.703 0.75 
2.48 - - - - 
2.715 7.36 8.19 Td Lod 
II - - - - 
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Leather 
10% 5% 
2.82 2.632 


10% 5% 
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Nylon fiber cuttings being light in weight were 
mixed in the proportion of 0.3%. The design mix 
for M25 concrete was determined as per the 
specifications of IS codes (IS:456, 2000; IS:10262, 
1982). The details of various design mixes are 
reported in Table 3. 


2. Soil and Solid Waste 

Subgrade soil was procured from 1 m depth below 
ground level and foundation soil was collected 
beyond | m depth. CKD was collected from Jaypee 
Cement Plant, Solan, Himachal Pradesh Cement 
kiln dust (CKD) contains some chemicals which are 
injurious to human health. Therefore, exposure of 
CKD to groundwater should be avoided. The 
groundwater table in Kurukshetra is at 25 m depth 
below ground level. Hence CKD was selected as an 
additive to stabilize subgrade soil only. Rice husk 
ash was procured from Kohinoor food limited, G.T 
Road, Murthal, Haryana, and Fly ash was collected 
from Thermal power plant, Assandh_ village, 
Panipat, Haryana. Rice husk ash (RHA) and Fly ash 
were selected to stabilize foundation soils. The 
physical properties of materials are reported in 
Table 4. All the tests were conducted on statically 
compacted soil samples (Jain and Puri, 2013b). 


ANALYSIS OF TEST RESULTS 

I. Concrete 

a) Weight of concrete: A significant decrease in the 
weight of concrete was observed when natural 
coarse aggregates were replaced with recycled 
aggregate. For plain cement concrete (PCC) weight 
was observed to be 7.92 kg. However, a lower range 
of weights was observed for concrete in which the 
natural aggregates were replaced with waste 
material. For concrete having recycled aggregates 
(RA) as constituents, weight ranges from 7.53 kg to 
7.78 kg. This can be attributed to the lower specific 
gravity of recycled aggregates. Also, a decreasing 
trend was observed for concrete in which fine 
aggregates were replaced with PVC waste, iron 
waste, and leather wastes. In the case of PVC 
concrete weight ranges from 7.65 kg to 7.87 kg and 
for leather concrete it ranges from 6.27 kg to 7.08 
kg. Hence lightweight concrete can be designed by 
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replacing natural aggregates with suitable waste 
materials. 


b) Slump: A decreasing trend was observed in the 
values of slump when natural aggregates were 
replaced by waste materials. For plain cement, the 
concrete value of the slump was observed to be 94 
mm. Values of slump decrease from 86 mm to 78 
mm for concrete in which coarse aggregates were 
replaced with recycled aggregate. For leather 
concrete slump ranges from 79 mm to 82 mm. A 
decrease in weight was observed with an increase in 
the percentage of pulverized leather. However, in 
the case of PVC concrete slump increases from 91 
mm to 92 mm with an increase in the percentage of 
PVC. Similarly, a decreasing trend in the values of 
the slump was observed for concrete having iron 
waste and nylon fibre cuttings. 


Observed values ofthe slump for a 
combination of 20% iron waste and 03% nylon fibre 
cuttings, 20% iron waste alone and 0.3% nylon fibre 
cuttings alone were observed as 75mm, 88mm and 
80 mm respectively. Better workability was 
observed in the case of PVC concrete as compared 
to the concrete having other wastes aggregates as 
constituents. This decreasing trend in workability 
can be attributed to the higher absorption capacity of 
waste materials. 


c) Compressive strength: Several cubes (a = 150 
mm) were prepared and were tested for compressive 
strength as per the specifications of IS codes [8, 9]. 
The average 28 days compressive strength of M25 
PCC was observed to be 28.72 kN/m*. When 
conventional aggregates were replaced with 
recycled aggregates in proportions of 100% to 25% 
the average compressive strength increased from 
31.28 kN/m? to 33.14 kN/m?. The achieved higher 
strength can be attributed to better bonding between 
recycled aggregates and conventional aggregates. 


When fine aggregates were replaced PVC 
waste in proportions of 10% to 5% the average 
compressive strength varies from 26.76 kN/m/? to 
28.87 kN/m*. 
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Table3: Details of Design Mix Calculations 


Mix Design 
Type of concrete 
Water (Litre) Cement (kg) Sand (kg) Coarse Aggregates (kg) 

PCC 0.46 1 1.34 2.86 
100% Recycled Aggregates 0.59 1 1.34 2.48 
50% Recycled Aggregates 0.58 1 1.34 2.67 
25% Recycled Aggregates 0.58 1 1.34 2:75 
5% PVC 0.46 1 1.34 2.86 
10% PVC 0.46 1 1.34 2.86 
5% Leather Waste 0.46 1 1.34 2.86 
10% Leather Waste 0.46 1 1.34 2.86 
20% Iron waste 0.53 1 1.18 222 
20% Iron waste + 0.3% Nylon fiber cuttings 0.53 1 1.18 222 
0.3% Nylon fiber cuttings 0.53 1 1.18 222 


Table 4: Physical Properties of Soil and Materials 


Materials 

Physical Properties 

Foundation Soil Subgrade Soil | Cement kiln dust Rice husk ash Fly ash 
Gravel (%) 0 - - 0 0 
Sand (%) 6.75 - 100 Zone i ie) 
Clay + silt (%) 93.25 100 - 76.78 O25 
Specific gravity 2.48 2.36 Doe, 1.95 2.09 
Shrinkage Limit 36.52 - - - - 
Liquid limit 54 45 NP 
Plastic limit 25 22 NP NP 
Plasticity index 29 22 
Is classification CH CI SM ML ML 
OMC (%) 14.71 18 - - - 
MDD (g/cc) 1.93 1.69 - - - 
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However, PVC concrete has shown lower strength 
in comparison with PCC but the characteristic 
strength was still achieved. This can be attributed to 
poor bonding between fine aggregate and PVC 
waste. 


When fine aggregates were replaced 
leather waste in proportions of 10% to 5% the 
average compressive strength varies from 0.938 
kN/m? to 9.07 kN/m?. Less than one-third of 
characteristic strength was achieved in this case. 
This can be attributed to the negligible bonding 
between leather waste and fine aggregates. 


Higher strengths were observed when 20% 
iron waste and 0.3% nylon wastes were used alone 
as well as in combination with other ingredients in 
concrete. In this case, 28 days strength varies from 
34.21 kN/m? to 37.92 kN/m’. This can be attributed 
to higher individual strengths of iron waste and 
nylon fibre cuttings. The variation of compressive 
strength with the age of concrete is shown in Figure 
1. 


2. Soil 

a) Shrinkage Characteristics: Shrinkage limit tests 
were conducted on foundation soil samples 
stabilized with Fly ash and Rice husk ash (RHA) 
(IS:2720-Part 6, 1972). The effect of these 
municipal solid wastes on the shrinkage of 
foundation soil was observed to be insignificant. 
However, a general increase in shrinkage limit from 
36.52 % to 43 % was observed as the percentage of 
RHA increased from 5% to 25%. However, a small 
decrease in the value of shrinkage limit from 36.52 
% to 35.65 % was observed for soil stabilized with 
fly ash. Reduction in shrinkage can be attributed to 
the reactions between soil and free lime content of 
RHA. Variation between shrinkage limit and 
percentage of additive is shown in Figure 2. 


b) Compaction Characteristics: Standard Proctor 
tests have been conducted to determine optimum 
moisture content (OMC) and maximum dry density 
(MDD) of both soils stabilized with various 
percentages of Fly ash, Rice husk ash (RHA), and 
Cement kiln dust (CKD) (IS:2720-Part 7, 1974). 
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Values of MDD and OMC for foundation soil have 
been observed as 1.93 g/cc and 14.91 % 
respectively. For subgrade soil, MDD and OMC 
have been observed as 1.69 g/cc and 18 %. It has 
been observed that MDD increases from 1.682 g/cc 
to 1.768 g/cc and OMC decreases from 19 % to 15 
% when subgrade soil was stabilized with CKD. 
This can be attributed to the cement content present 
in CKD. However, an opposite trend has been 
observed for foundation soil stabilized with Fly ash 
and RHA. For Fly ash, stabilized samples, MDD 
decreases from 1.8 g/cc to 1.77 g/cc as the 
percentage of Fly ash increases. But OMC increases 
from 12.9 % to 22.5 % with an increase in the 
percentage of Fly ash. Similarly, for RHA stabilized 
samples, MDD decreases from 1.81 g/cc to 1.54 g/cc 
as the percentage of RHA increases. But OMC 
increases from 19.4 % to 30.23 % with an increase 
in the percentage of RHA. This can be attributed to 
the lower specific gravity and high absorption 
capacity of Fly ash and RHA. The variation of MDD 
and OMC with the percentage of additive is shown 
in Figure 3 and Figure 4. 


c) Strength Parameters: Unconfined compressive 
strength tests were performed in order to determine 
strength characteristics of soils stabilized with 
various municipal solid wastes (IS:2720-Part 10, 
1973). Values for shear strength of foundation soil 
and subgrade soil have been observed as 3.357 
N/cm? and 0.85 N/cm? respectively. It has been 
observed that the shear strength of subgrade soil 
stabilized with CKD increases from 0.85 N/cm/? to 
2.9 N/cm?* with an increase in the percentage of 
CKD. Shear strength of foundation soil stabilized 
with Fly ash increases from 3.36 N/cm* to 3.65 
N/cm/ with an increase in the percentage of Fly ash. 
Similarly, the Shear strength of foundation soil 
stabilized with RHA increases from 3.36 N/cm? to 
3.85 N/cm? with an increase in the percentage of 
RHA. This increase in the shear strength of soil can 
be attributed to the reactions between soil and free 
lime content of these solid wastes. Significant 
improvements in other engineering properties have 
also been reported by various authors (Jain and Puri, 
2013a). Variation of UCS with the percentage of 
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additive for different soil samples is shown in Figure 


5. 


CONCLUSIONS 


Based on the results, following conclusions have 


been drawn: 


1. 


Recycled aggregates are better alternatives to 
conventional coarse aggregates. A significant 
increase in the compressive strength of concrete 
has been observed in the case of recycled 
aggregates, iron waste, and nylon fiber cuttings. 
However, a lower slump has been observed in 
these cases. It is recommended to use 
admixtures for improving the workability of 
resultant concrete. 

Municipal solid wastes like PVC waste, nylon 
fiber cuttings, and steel scrap can only be used 
in small quantities because of their lower 
specific gravities and incompatibility with 
earthen materials. Also, resistance against fire 
is not assured for concrete having these wastes 
as constituents. 

Leather waste should not be used as a 
replacement for conventional aggregates. 

It is recommended by authors to use CKD as a 
stabilizer for problematic soil at shallow depths 
and rice husk ash for soils beyond shallow 
depths. 

These industrial wastes do consist of harmful 
chemicals, their exposure to the human 
atmosphere and groundwater table must be 
avoided. 
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ABSTRACT 


Landfilling is the most common technique used in municipal solid waste (MSW) management, In India, low 
laying areas in the vicinity of urban expanses are preferred for landfills. A site in Warriana area of Jalandhar 
district of Punjab (India) has been selected for the present study. Jalandhar is facing a big problem with solid 
waste management and its disposal. A total of 490-550 TPD of MSW is being generated in Jalandhar and dumped 
at Warriana Site. The waste when comes in contact with water produces a hazardous substance called “Leachate’’. 
This leachate, if not isolated properly, penetrates into the ground and pollutes the groundwater. In this study, the 
pollution potential of the leachate present at the Warriana site has been assessed. The possible utilization of waste 
tire scraps and aggregates mixed in a 1:1 ratio for leachate re-circulation has been explored. A reduced pollution 


potential of the leachate has been observed. 


Keywords: Leachate, solid waste, tire scraps, leachate recirculation, landfill 


INTRODUCTION 


Today, economical management of the MSW and 
the leachate generated from it has become a global 
challenge (USEPA 1998). Increasing population, 
rapid urbanization and worldwide industrialization 
are the major contributing factors. In the State of 
Punjab, most of the cities are dealing with this 
problem. This can be resolved by the use of various 
waste management alternatives such as land filling. 
However, authorities are facing huge difficulties to 
manage such an enormous amount of waste because 
of incorrect design, improper guidelines and lack of 
research in this area. 


In Punjab, nearly 4500 tons of municipal 
solid waste is generated daily. Out of this, the 4 
municipal corporations (Ludhiana, Jalandhar, 
Amritsar, and Patiala) account for the generation of 
2000 Tons per day (TPD) of MSW in Punjab 
(Kaushik et al. 2014). Out of this, only Jalandhar 


@ Corresponding Author’s Affiliation: 
Department of Civil Engineering 


Lyallpur Khalsa College of Engineering, Jalandhar, Punjab (India) 


city generates about 550 TPD of MSW (Kaushik and 
Sethi, 2007). 


Similarly, scrap tires are also another waste 
material which require proper collection and the 
disposal. Automobile tires are generally made of 
natural rubber elastomers, polymers and _ other 
additives. Although their composition varies with 
the type of the tire and the manufacturer company. 
Steel reinforcement are also provided in the tires to 
improve its strengths. These tires even after use, 
maintain their chemical composition. They require 
thousands of years to fully decompose. It has been 
estimated that over 200 million scrap tires are 
generated in India on annual basis. Only 18% of 
them are safely reused or recycled. Even through 
tires comprises only 1% of the total waste generated 
in India, they present significant problems due to 
their non-biodegradable nature. Even if we set tires 
on fire, they produce a lot of air pollution releasing 
significant amount of pollutant gases in the air. 
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The use of tire scraps and aggregates would 
allow the leachate to drain/flow from the soil, so that 
it is treated and does not pollute the underlying soil 
and groundwater. Also, the utilization of scrap tires, 
would resolve the problem of their disposal. The 
scrap tire later purifies the leachate in the landfill 
and also prevents buildup of the hydraulic head 
within the filter. The objective of the study was to 
check the performance of the scrap tire-gravel 
mixture in treatment systems of leachate in landfills 
and to highlight their usefulness as a drainage 
medium. 


STUDY AREA 


The study area Jalandhar lies in Doaba region of 
Punjab. It is located between two rivers. The whole 
state is divided into 7 clusters. Jalandhar cluster 
contains 26 municipal committees including 
Municipal Corporation Jalandhar. Solid waste 
management in Jalandhar is under the Sanitation 
Department and health section department of the 
Municipal Corporation Jalandhar. 


MATERIALS AND METHODOLOGY 


Materials used in the study are as follows: 


1. Tire Scraps: The tire chips used for this study 
were obtained from the material suppliers in the 
local market of Jalandhar city. Conventional 
grain size distribution analysis (ASTM: D422- 
63, 1998) was performed on the tire chips used 
in the study, due to their irregular sizes. The 
sizes of tire chips ranged from 25-50 mm in 
length, with chip size of an average of 40 mm. 
The test specimen comprised of wire free tire 
chips. The specific gravity (Gs) of tire shreds 
ranges from 1.032 to 1.36, depending upon 
amount of steel wires in the tire (Kaushik et al., 
2010). Reddy and Saichek (1998) found that the 
hydraulic conductivity of tire shreds is reduced 
to just 0.01 cm/second under a normal stress of 
21000 psf. If the tire shreds are compressed by 
65% under this level of normal stress, they met 
the minimum hydraulic — conductivity 

requirement for drainage materials which is 

0.001cm/s. 
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Fig: 1 Tire Scraps 


2. The aggregates used in this study were recycled 
concrete aggregates (RCA). They are obtained 
by crushing the waste concrete cubes at the 
concrete laboratory of the institute. The size 
chosen for this study was 40 mm to 20 mm. 





Fig: 2 Recycled Aggregates 


Experimental Setup: In order to create conditions 
like the landfill, a 12-inch pipe was taken with its 
one end closed and the other end open. A tap was 
also attached to the bottom of the pipe for the 
extraction of the leachate. The RCA and tire scraps 
were laid in_ layers, keeping each layer 
approximately 4 to 5 inches thick. The leachate 
which has been obtained from the Warriana site is 
recirculated several times and its properties are 
checked after every 10 cycles for upto 30 cycles. 
The various parameters for which the leachate was 
tested were: 


e BOD (Bio-Chemical Oxygen Demand) 
e COD (Chemical Oxygen Demand) 

e Chlorides 

e Hardness 

e Acidity 

e Total Solids 

e Total Suspended Solids 
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RESULTS AND DISCUSSION 


The leachate sample was taken from the Warriana 
site. Then it was tested initially and after particular 
no. of cycles. During this, the reduction in the BOD 
and the COD values of the leachate was observed. 
The improvement observed in the quality of the 
leachate from its initial stage to up to 30 cycles from 


the filter media have been presented in Figure 1 to 
A. Chloride (in mg/l) 


B initial leachate & ist cycle 


@ 2nd cycle 3rd cycle 





Fig.3: Chlorides observed till 30 Cycles 


BOD (in mg/l) 


= Initial leachate = 1st cycle 
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Fig.l: BOD observed till 30 Cycles 
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Fig.4: Acidity observed till 30 Cycles 
COD (in mg/l) 


B initial leachate = ist cycle 
= 2nd cycle = 3rd cycle 





Fig.2: COD observed till 30 Cycles 
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Parameter 

Biochemical oxygen 
demand (BOD) 

Chemical oxygen demand 


(COD) 


Chloride 


Acidity 


Hardness 


Total solids 


Total suspended solids 


Table 1: Change in pollution potential of leachate after passing through 30 cycles 


Results of initial leachate 


180 mg/l 


11.2 mg/l 


464.85 mg/l 


66 mg/l 


940 mg/l 


292 mg/l 


17 mg/l 
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Results after 30 cycles 


115 mg/l 


3.2 mg/l 


204.93 mg/l 


18 mg/l 


610 mg/l 


10 mg/l 


1.0 mg/l 


Observations 


2.77% decrease in the value of BOD 


26.78 % decrease in the value of COD 


36.56 % decrease in the value of 
chloride 


66.66 % decrease in the value of 
acidity 


17.02 % decrease in the value of 
hardness 


68.15 % decrease in the values of total 
solids 


70.58 % decrease in the values of total 
suspended solids 
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CONCLUSION 


The leachate collected from Warriana dumping site 
of Jalandhar has been characterized for different 
parameters such as Biochemical Oxygen Demand 
(BOD), Chemical Oxygen Demand (COD), 
Hardness, Chlorides, Sulfates, Total solids and 
dissolved solids to check its potential to pollute the 
underground aquifers. It was analyzed that leachate 
from this site is of moderate strength and requires 
some treatment before its safe disposal either on 
land surface or in water bodies. For this the leachate 
has been passed through the layers of gravel and tire 
chips. 


After passing the leachate by 30 cycles the 
pollution potential measured as BOD is reduced by 
36.11%, COD is reduced by 72%, Chloride is 
reduced by 56%, acidity is reduced by 73%, 
Hardness is reduced by 36%, total solids are reduced 
by 97 % and total suspended solids are reduced by 
95%. The analysis shows that the tire scraps and 
RCA can be beneficially reused to reduce the 
pollution problems arising from the landfill. The 
approach used in this study can be used for other 
landfill sites to reduce pollution potential of the 
leachate. 
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ABSTRACT 


A preliminary assessment of earthquake hazards for the district of Pathankot (Punjab) has been carried out. The 
study area has already experienced strong ground motions of intensity VII in the past during the 1905 Kangra 
earthquake. The epicenter of the 1905 Kangra earthquake was only 60 km away from the study area. Hence, the 
proximity of the study area to the northwestern Himalayas makes it susceptible to damage due to large earthquakes 
in the future. Therefore, deterministic seismic hazard analysis (DSHA) has been carried out for the district. For 
this purpose, earthquake data has been collected from various seismological agencies e.g. NDMA, IMD, ISC- 
UK, and USGS, and a comprehensive earthquake catalog has been prepared. Several ground motion prediction 
equations (GMPEs) have been reviewed and analyzed to select a GMPE appropriate for carrying out DSHA. 
Earthquake hazard parameters have been formulated in terms of peak ground accelerations (PGA). Based on the 
results, a seismic hazard map for the district has been prepared. The study shows that the region can experience 
severe strong ground motions originating from northwestern Himalayan earthquakes and hence, structures must 
be designed based on hazard parameters determined considering the local tectonic setting. The developed hazard 
map would help engineers and town planners in designing earthquake-resistant structures and retrofitting existing 
Structures. 


Keywords: DSHA, Maximum magnitude potential, Peak ground acceleration, Pathankot 


INTRODUCTION for the government to take suitable mitigation 


measures. 
Earthquakes are the most hazardous natural events 


which can inflict massive damage to structures and 
can wipe out entire cities in minutes. In the past few 
decades, the earthquakes in India have gone up as 
the Indian plate and Eurasian plate are bashing at a 
rate of 50 mm/year (Peltzer and Saucier, 1996). The 
mighty Himalayas are born due to the clashing of 
these two plates, which makes the terrain area in the 
vicinity susceptible to moderate to severe 
earthquakes. Some of the significant seismic events 
that occurred in this area are the Uttarkashi 
earthquake (1991), the Chamoli earthquake (1999), 
and the Nepal earthquake (2015). The wreckage 
caused by these earthquakes has been a wake-up call 
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Seismic zoning map of India as per BIS- 
1893: Part-1 (2016) offers a crude idea of the 
earthquake hazard as it gets updated only after the 
occurrence of an earthquake. Hence, for an 
emerging country like India, appropriate steps for 
seismic hazard assessment are necessary to estimate 
an optimum and reliable value of possible 
earthquake ground motion. These prognosticated 
values could be used as an input to assess the seismic 
vulnerability of an area, based on which new 
construction and the restoration works of existing 
structures can be carried out. The purpose of a 
seismic hazard analysis (SHA) is to quantify 
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potential damage and loss due to future earthquakes. 
It includes a quantitative evaluation of ground 
shaking hazards at a particular site. SHA can be used 
to make regional macro or micro zoning maps which 
are useful in earthquake-resistant building design 
and construction, land use planning, emergency 
plan, and estimation of possible economic loss. 


Two perspectives, probabilistic seismic 
hazard analysis (PSHA) and deterministic seismic 
hazard analysis (DSHA) are commonly used for 
seismic hazard assessment. In DSHA, a specific 
earthquake scenario is assumed based on both 
historical and instrumental data, and hazard 
parameters are estimated considering the attenuation 
characteristics of the region. In this method, the 
source causing the maximum hazard at any location 
is considered to be the controlling source that poses 
a maximum threat to that particular location. 
Probabilistic seismic hazard analysis (PSHA) is a 
methodology that considers the uncertainties in size, 
location, and occurrence of an earthquake, and 
estimates the likelihood of various levels of 
earthquake-caused ground motions being exceeded 
in a given time period at a location. The results of 
such an analysis are expressed as_ estimated 
probabilities per year or estimated annual 
frequencies. 


The results of the seismic hazard analysis 
are formulated generally in terms of peak ground 
acceleration (PGA) and spectral acceleration (Sa). 
Conventionally, PGA is used to quantify ground 
motion and used to calculate lateral forces and shear 
stresses in the equivalent-static-force procedures in 
building codes, and in liquefaction analysis. The Sa 
is the maximum acceleration experienced by a 
damped single-degree-of-freedom oscillator and 1s a 
crude representation of building response. These 
parameters are the prime inputs to seismic hazard 
analysis and hence must be calculated with extreme 
caution. 


In the present study, an attempt has been 
made to assess seismic hazards for the District of 
Pathankot in State of Punjab using a deterministic 
approach. The results of the study would help in the 
more efficient design of earthquake-resistant 
structures and the planning of rescue arrangements. 
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STUDY AREA 


Pathankot district is located in the north zone of the 
State of Punjab (India). The district is spread in an 
area of 929 km7?. It is located at the foothills of 
Shivalik Hills. It shares international borders with 
the Narowal District of Pakistani Punjab and state 
borders with the Kathua District of Jammu & 
Kashmir and Chamba and Kangra districts of 
Himachal Pradesh. The two main rivers, Beas and 
Ravi, passes through the district. The Latitude is 
32.324276° N and the longitude is 75.59063° E. 


The District is only 60 km away from the 
epicenter of the 1905 Kangra earthquake. The 1905 
Kangra earthquake occurred in the Kangra Valley 
and the Kangra region of the Punjab Province 
(modern-day Himachal Pradesh) in India on 4" 
April 1905. The earthquake measured 7.8 on the 
surface wave magnitude scale and killed more than 
20,000 people. Some of the major earthquakes that 
occurred near the region during the last 200 years 
are the 1934 Bihar-Nepal earthquake (Mw 8.4), 
1950 Assam earthquake (Mw 8.7), 1905 Kangra 
earthquake (Mw 7.8), 1991 Uttarkashi (Mw 6.8), 
1993 Killari earthquake (Mw 6.2), 1999 Chamoli 
earthquake (Mw 6.8), 2005 Kashmir earthquake 
(Mw 7.6) and 2015 Gorkha earthquake (Mw 7.8). 


An area in 300 km radius around the 
Pathankot district has been selected as the seismic 
study region (SSR). 


EARTHQUAKE HISTORY OF THE REGION 


A comprehensive earthquake catalogue is a 
prerequisite for earthquake hazard analysis. A 
reliable seismic hazard assessment for a region 
strongly depends on the data statistics of the events. 
There are two types of database available: pre- 
instrumental and instrumental period database. 
Records of pre-instrumental events are of immense 
importance for the compilation of earthquake 
catalogue. Pre-instrumental data are available for 
damaging earthquakes and have been collected from 
catalogue of National Disaster Management 
Authority of India (NDMA, 2011). The instrumental 
data have been collected from various national and 
international seismological agencies such as 
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National Disaster Management Authority (NDMA), 
Indian Meteorological Department (IMD), 
International Seismological Centre (ISC-UK) and 
United States Geological Survey (USGS). A 
comprehensive earthquake catalogue has been 
prepared for a period from 1800 AD to 2017 AD by 
combining the pre-instrumental and instrumental 
data collected. The compiled catalogue has the 1256 
events of magnitude Mw > 4 occurred in the seismic 
study region (between 28° - 33° N to 73° - 77° E) till 
December 2017. An Epicenter map has been 
developed using the prepared catalogue as shown in 
Figure 1. 


4-4.9 Mw 
5-5.9 Mw 
6-6.9 Mw 
7-7.9 Mw 


*erClle 


Pathankot District 





(N) 


Latitude 














Longitude (E) 


Fig. 1: Epicentral Map of Study Area 


The events in the catalogue have been 
converted to a common scale of moment magnitude 
(Mw). For this, numerous empirical relations are 
available to convert different magnitude scales to 
moment magnitude (Mw). Pre-instrumental data is 
usually available on MMI scale and have been 
converted to moment magnitude (Mw) using the 
Gutenberg-Richter (1956) relation. The following 
correlations developed by Scordilis (2006) between 
mb - Mw and Ms - Mw have been used. For the 
conversion of local magnitude (ML) to moment 
magnitude (Mw), regional correlation is generally 
preferred. Hence, the following correlation 
developed by Kolathayar et al. (2012) for India and 
adjoining areas has been used. The following 
correlation developed by Yenier et al. (2008) 
between Md - Mw has been used. 
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SEISMOTECTONICS OF THE REGION 


A tectonic map developed for the seismic study 
region by Puri and Jain (2018) has been used in this 
study (see Figure 2). The tectonic map was 
developed using Seismotectonic Atlas of India and 
its Environs (SEISAT) (Dasgupta et al., 2000). In 
order to identify the active tectonic sources likely to 
generate strong ground motions in the study area, 
the earthquake data from the catalogue have been 
superimposed on to the tectonic map. A total of 10 
tectonic features which are associated with 
earthquakes of moment magnitude 4.0 and above 
have been identified as active tectonic sources, out 
which 08 sources have been selected as potential 
seismogenic sources. The values of maximum 
observed magnitude (Mobs) for various 
seismogenic sources have been reported in Table 1. 


Table 1: Maximum Observed Magnitudes for Various Sources 


S.No. Seismogenic Source Mobs 
1 Sargodha Lahore Delhi Ridge (SLDR) 6.5 
2 Main Boundary Thrust (MBT) 8 
3 Main Frontal Thrust (MFT) 5 
4 Jwala Mukhi Thrust (JMT) J0 
5 Main Crustral Thrust (MCT) Tia 
6 Sundar Nagar Fault (SNF) | 
7 Kaurik Fault System (KFS) 6.8 
8 Lineament System of SLDR (LSLDR) 6 


It has been observed that the district of 
Pathankot has always been under the threat of big 
earthquakes due to its proximity with Main 
Boundary Thrust (MBT) and Main Frontal Thrust 
(MFT) of Himalayan Thrust System. 


ESTIMATION OF MAXIMUM MAGNITUDE 
POTENTIAL FOR THE STUDY REGION 


Damaging earthquakes of magnitude (Mw) of 5 and 
above have been considered for estimation of 
maximum magnitude potential (Mmax). The Mumax 
values for the potential se1smogenic sources present 
in the study area have been estimated by adding an 
increment 0.5 to the respective Mops values (Gupta, 
2002). The Mmax values in terms of moment 
magnitude estimated for the seismic study area 
range from 5.5 to 8.5.The Mmax and total fault 
lengths (TFL) values for the selected seismogenic 
sources have been reported Table 2. 
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Fig. 2: Tectonic setup of the seismic study region 


Table 2: Maximum Magnitude Potential for Various Sources 


Seismogenic Total Fault Maz 
Source Length (TFL) in 
kms 
SLDR Area Source 7 
MBT 825 8.5 
MFT 32 5.5 
JMT 290 6.0 
MCT 769 7.8 
SNF 101 75 
KFS its i 7.3 
LSDSR 97 6.5 


ESTIMATION OF SEISMIC HAZARD 


In this study, a global ground motion prediction 
equation (GMPE) for shallow crustal earthquakes 
developed by Graizer and Kalkan (2007) has been 
used for estimating seismic hazard for the district of 
Pathankot. For the computation of hazard at a grid 
point, ground motion parameters have been 
computed w.r.t various seismogenic sources and the 
source causing maximum ground motion at the point 
of interest has been identified. The maximum 
magnitude potential of that seismogenic source has 
been considered as the controlling earthquake for 
the grid point under consideration. The maximum 
PGA value, i.e. PGAmax, for the controlling 
earthquake at a grid point has been considered to 
prepare a seismic hazard map of the district. The 
value of PGAmax for all grid points ranges from 
0.254¢ to 1.89g. A deterministic seismic hazard map 
has been developed for the district using nearest 
neighbor interpolation technique and is shown in 
Figure 3. 


Annals of Earth, Environment, and Infrastructure Engineering 


Volume 1, Issue 1 (November-December), pp. 33-37 





32.6 


Latitude (N) 
Ww 
No 
Ww 
T 








52 75.4 . 75.6 75.8 76 
Longitude (E) 


Fig. 3: Deterministic Seismic Hazard Map for Pathankot 


A very high to severe seismic hazard has 
been observed for the entire district. The expected 
PGA for north eastern part of Pathankot is as high as 
1.89g. It has been observed that the PGA reported 
for this region in IS:1893-Part 1 (2016) quite 
underestimates the expected seismic hazard. 


LIMITATION OF THE STUDY 
The study has following limitations: 


1. Seismic hazard maps based on probabilistic 
approach would give more realistic seismic 
hazard scenario. 

2. Due to preliminary nature of the study, the Mmax 
values have been calculated using incremental 
method only. However, a better assessment can 
be made using some advanced methods, e.g., 
Mark (1997), Wells and Coopersmith (1994) 
and Kijko (2011). 

3. The catalogue should be screened _ for 
foreshocks and aftershocks (Gardner and 
Knopoff, 1974) and completion periods of 
should be determined using CUVI method 
(Mulargia and Tinti, 1985) and Stepp’s (Stepp 
1972) methods. 

4. The study been carried out for rock sites and 
hence it is necessary to carry out ground 
response analysis (e.g. Puri and Jain, 2021) for 
the sites underlain by soils in order to account 
for the nonlinearity of soils at severe peak 
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ground accelerations. 

5. A detailed study is required for this region of 
Punjab so that the hazard parameters given in 
IS:1893-Part 1 (2016) could be suitably revised. 


CONCLUSIONS 


In the present study, the seismic hazard analysis 
using deterministic approach has been carried out 
for estimating peak ground acceleration (PGA) at 
rock sites for district of Pathankot. It has been 
observed that the district is prone to high to severe 
earthquake ground motions ranging up to 1.89g. 
Hence, a comprehensive evaluation of earthquake 
hazard for this region is must. The results of this 
study would help and motivate engineers and 
investigators in earthquake resistance design of 
structures and planning of mitigation measures. 
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